skip to main content


Search for: All records

Creators/Authors contains: "Zrake, Jonathan"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Abstract

    We show that gas disks around the components of an orbiting binary system (so-called minidisks) may be susceptible to a resonant instability that causes the minidisks to become significantly eccentric. Eccentricity is injected by, and also induces, regular impacts between the minidisks at roughly the orbital period of the binary. Such eccentric minidisks are seen in vertically integrated, two-dimensional simulations of a circular, equal-mass binary accreting from a circumbinary gas disk with a Γ-law equation of state. Minidisk eccentricity is suppressed by the use of an isothermal equation of state. However, the instability still operates and can be revealed in a minimal disk-binary simulation by removing the circumbinary disk and feeding the minidisks from the component positions. Minidisk eccentricity is also suppressed when the gravitational softening length is large (≳4% of the binary semimajor axis), suggesting that its absence could be an artifact of widely adopted numerical approximations; a follow-up study in three dimensions with well-resolved, geometrically thin minidisks (aspect ratios ≲0.02) may be needed to assess whether eccentric minidisks can occur in real astrophysical environments. If they can, the electromagnetic signature may be important for discriminating between binary and single black hole scenarios for quasiperiodic oscillations in active galactic nuclei; in turn, this might aid in targeted searches with pulsar timing arrays for individual supermassive black hole binary sources of low-frequency gravitational waves.

     
    more » « less
  2. Abstract

    Observations and theory suggest that core-collapse supernovae can span a range of explosion energies, and when sub-energetic the shockwave initiating the explosion can decelerate to speeds comparable to the escape speed of the progenitor. In these cases, gravity will complicate the explosion hydrodynamics and conceivably cause the shock to stall at large radii within the progenitor star. To understand these unique properties of weak explosions, we develop a perturbative approach for modeling the propagation of an initially strong shock into a time-steady, infalling medium in the gravitational field of a compact object. This method writes the shock position and the post-shock velocity, density, and pressure as series solutions in the (time-dependent) ratio of the freefall speed to the shock speed, and predicts that the shock stalls within the progenitor if the explosion energy is below a critical value. We show that our model agrees very well with hydrodynamic simulations, and accurately predicts (for example) the time-dependent shock position and velocity and the radius at which the shock stalls. Our results have implications for black hole formation and the newly detected class of fast X-ray transients (FXTs). In particular, we propose that a “phantom shock breakout”—where the outer edge of the star falls through a stalled shock—can yield a burst of X-rays without a subsequent optical/UV signature, similar to FXTs. This model predicts the rise time of the X-ray burst,td, and the mean photon energy,kT, are anticorrelated, approximately asTtd5/8.

     
    more » « less
  3. ABSTRACT

    The upcoming Laser Interferometer Space Antenna (LISA) is expected to detect gravitational waves (GWs) from massive black hole binaries (MBHB). Finding the electromagnetic (EM) counterparts for these GW events will be crucial for understanding how and where MBHBs merge, measuring their redshifts, constraining the Hubble constant and the graviton mass, and for other novel science applications. However, due to poor GW sky localization, multiwavelength, time-dependent EM models are needed to identify the right host galaxy. We studied merging MBHBs embedded in a circumbinary disc (CBD) using high-resolution two-dimensional simulations, with a Γ-law equation of state, incorporating viscous heating, shock heating, and radiative cooling. We simulate the binary from large separation until after merger, allowing us to model the decoupling of the binary from the CBD. We compute the EM signatures and identify distinct features before, during, and after the merger. Our main result is a multiband EM signature: we find that the MBHB produces strong thermal X-ray emission until 1–2 d prior to the merger. However, as the binary decouples from the CBD, the X-ray-bright minidiscs rapidly shrink in size, become disrupted, and the accretion rate drops precipitously. As a result, the thermal X-ray luminosity drops by orders of magnitude, and the source remains X-ray dark for several days, regardless of any post-merger effects such as GW recoil or mass-loss. Looking for the abrupt spectral change where the thermal X-ray disappears is a tell-tale EM signature of LISA mergers that does not require extensive pre-merger monitoring.

     
    more » « less
  4. Abstract Linear analysis of gas flows around orbiting binaries suggests that a centrifugal barrier ought to clear a low-density cavity around the binary and inhibit mass transfer onto it. Modern hydrodynamics simulations have confirmed the low-density cavity, but show that any mass flowing from large scales into the circumbinary disk is eventually transferred onto the binary components. Even though many numerical studies confirm this picture, it is still not understood precisely how gas parcels overcome the centrifugal barrier and ultimately accrete. We present a detailed analysis of the binary accretion process, using an accurate prescription for evolving grid-based hydrodynamics with Lagrangian tracer particles that track the trajectories of individual gas parcels. We find that binary accretion can be described in four phases: (1) gas is viscously transported through the circumbinary disk up to the centrifugal barrier at the cavity wall, (2) the cavity wall is tidally distorted into accretion streams consisting of near-ballistic gas parcels on eccentric orbits, (3) the portion of each stream moving inwards of an accretion horizon radius r ¯ ≃ a —the radius beyond which no material is returned to the cavity wall—becomes bound to a minidisk orbiting an individual binary component, and (4) the minidisk gas accretes onto the binary component through the combined effect of viscous and tidal stresses. 
    more » « less
  5. Abstract Core-collapse supernovae can display evidence of interaction with preexisting, circumstellar shells of material by rebrightening and forming spectral lines, and can even change types as hydrogen appears in previously hydrogen-poor spectra. However, a recently observed core-collapse supernova—SN 2019tsf—was found to brighten after roughly 100 days after it was first observed, suggesting that the supernova ejecta was interacting with surrounding material, but it lacked any observable emission lines and thereby challenged the standard supernova-interaction picture. We show through linear perturbation theory that delayed rebrightenings without the formation of spectral lines are generated as a consequence of the finite sound-crossing time of the postshock gas left in the wake of a supernova explosion. In particular, we demonstrate that sound waves—generated in the postshock flow as a consequence of the interaction between a shock and a density enhancement—traverse the shocked ejecta and impinge upon the shock from behind in a finite time, generating sudden changes in the shock properties in the absence of ambient density enhancements. We also show that a blast wave dominated by gas pressure and propagating in a wind-fed medium is unstable from the standpoint that small perturbations lead to the formation of reverse shocks within the postshock flow, implying that the gas within the inner regions of these blast waves should be highly turbulent. 
    more » « less
  6. Abstract

    The stellar cataclysms producing astronomical transients have long been modeled as either a point-like explosion or jet-like engine ignited at the center of a spherically symmetric star. However, many stars are observed, or are expected on theoretical grounds, not to be precisely spherically symmetric, but rather to have a slightly flattened geometry similar to that of an oblate spheroid. Here we present axisymmetric two-dimensional hydrodynamical simulations of the dynamics of point-like explosions initiated at the center of an aspherical massive star with a range of oblateness. We refer to these exploding aspherical stars as “ellipsars” in reference to the elliptical shape of the isodensity contours of their progenitors in the two-dimensional axisymmetric case. We find that ellipsars are capable of accelerating expanding rings of relativistic ejecta. which may lead to the production of astronomical transients including low-luminosity gamma-ray bursts, relativistic supernovae, and fast blue optical transients

     
    more » « less
  7. null (Ed.)
  8. ABSTRACT During the final stages of a compact object merger, if at least one of the binary components is a magnetized neutron star (NS), then its orbital motion substantially expands the NS’s open magnetic flux – and hence increases its wind luminosity – relative to that of an isolated pulsar. As the binary orbit shrinks due to gravitational radiation, the power and speed of this binary-induced inspiral wind may (depending on pair loading) secularly increase, leading to self-interaction and internal shocks in the outflow beyond the binary orbit. The magnetized forward shock can generate coherent radio emission via the synchrotron maser process, resulting in an observable radio precursor to binary NS merger. We perform 1D relativistic hydrodynamical simulations of shock interaction in the accelerating binary NS wind, assuming that the inspiral wind efficiently converts its Poynting flux into bulk kinetic energy prior to the shock radius. This is combined with the shock maser spectrum from particle-in-cell simulations, to generate synthetic radio light curves. The precursor burst with a fluence of ∼1 Jy·ms at ∼GHz frequencies lasts ∼1–500 ms following the merger for a source at ∼3 Gpc (Bd/1012 G)8/9, where Bd is the dipole field strength of the more strongly magnetized star. Given an outflow geometry concentrated along the binary equatorial plane, the signal may be preferentially observable for high-inclination systems, that is, those least likely to produce a detectable gamma-ray burst. 
    more » « less
  9. null (Ed.)
  10. null (Ed.)